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ABSTRACT 

In this, the third in a series of three papers concerning the SuperCOSMOS Sky Survey, 
we describe the astrometric properties of the database. We describe the algorithms 
employed in the derivation of the astrometric parameters of the data, and demonstrate 
their accuracies by comparison with external datasets using the first release of data, 
the South Galactic Cap survey. We show that the celestial co-ordinates, which are tied 
to the International Celestial Reference Frame via the Tycho-2 reference catalogue, are 
accurate to better than ± 0.2 arcsec at J,R ~ 19,18 rising to ± 0.3 arcsec at J,R ~ 22,21 
with positional dependent systematic effects from bright to faint magnitudes at the 
~ 0.1 arcsec level. The proper motion measurements are shown to be accurate to 
typically ± 10 mas yr" 1 at J,R ~ 19,18 rising to ± 50 mas yr _1 at J,R ~ 22,21 and 
are tied to zero using the extragalactic reference frame. We show that the zeropoint 
errors in the proper motions are < 1 mas yr _1 for R > 17 and are no larger than 
~ 10 mas yr _1 for R < 17 mas yr -1 . 

Key words: astronomical databases: miscellaneous - catalogues - surveys - astrom- 
etry - reference systems - stars: proper motions 



1 INTRODUCTION 

In Paper I of this series (Hambly et al. 2001a) we describe 
the SuperCOSMOS Sky Survey (hereafter SSS)R. This am- 
bitious project ultimately aims to digitise the entire sky atlas 
Schmidt plate collections in three colours (BRI), one colour 
(R) at two epochs. Paper II in this series (Hambly, Irwin & 
MacGillivray 2001b) describes the image detection, param- 
eterisation, classification and photometric calibration tech- 
niques for the survey. In this, the third paper of the series, we 
describe in detail the calibration of the astrometric parame- 
ters contained within the survey data. Paper I is intended as 
a user guide to the SSS while Papers II and III provide more 
technical details concerning the derivation and calibration of 
the object parameters. The first release of data consisted of 
~ 5000 square degrees of the southern sky at high Galactic 
latitude (|6| > 60°), and is known as the South Galactic Cap 
survey. Although these three papers make explicit reference 
to SGC data, all details are generally applicable to SSS data 
at Galactic latitudes \b\ > 30°. At low Galactic latitudes, 
particularly towards the Galactic centre, crowding degrades 
astrometric and photometric performance. 



database available online at http://www-wfau.roe.ac.uk/sss 



The photographic plate material used for the SSS sur- 
vey consist of sky-limited Schmidt photographic glass plate 
and film originals, or glass copies of glass originals, taken 
with the UK, ESO and Palomar Oschin Schmidt Telescopes 
- for more details see Paper I, Morgan et al. (1992) and refer- 
ences therein. Leaving aside such details as object detection 
and parameterisation (see Paper II) the solution to the prob- 
lem of assigning celestial co-ordinates to each object requires 
comparison with a reference catalogue and the application 
of a plate model. Sc hmidt plate astrometry has a long his- 
tory (see Section 2.1 for some references) involving reference 
catalogues of increasing density and precision along with in- 
creasingly sophisticated plate models - for a concise review, 
see Morrison, Smart & Taff (1998). These authors question 
the need for a conventional plate model at all and advocate 
a novel technique making use of an 'infinitely overlapping 
circles filter'. However, we note that they concede that, for 
large-scale survey work at least, a global linear plate model 
plus non-linear correction via an empirical plate 'mask' is 
approp riate . It is just such a method that we describe in 
Section 2.1. In Section^, we proceed to test the astrometric 
solutions against CCD drift scan data encompassing plate 
boundaries - ie. a test of the true external errors, both ran- 
dom and systematic, as a function of magnitude and plate 
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position. Section g also discusses our results in comparison 
to extragalactic objects denning the International Celestial 
Reference Frame (Ma et al. 1998). 

In addition to providing celestial co-ordinates, the SSS 
survey catalogues include proper motions. In the south- 
ern hemisphere, all objects paired between the SERC-J/EJ 
(first epoch, hereafter Bj) and SERC-ER/AAO-R (second 
epoch, hereafter R) pla tes have a proper motion measure- 
ment (see Section 2.2.1). The methodology for determining 



the proper motions closely follows that described in Eva ns fc 
Irwin (1995; hereafter EI95) and is detailed in Section 
Small changes made necessary by the plate material used 
are also described. Section ^| contains a comparison between 
our derived proper motion measurements against an exter- 
nal dataset, and we analyse in some detail the relative con- 
tribution of random and systematic errors on these proper 
motions as a function of magnitude. 

As we will describe in Section 2.1 the celestial co- 
ordinates are tied to the Hipparcos reference frame via the 
Tycho-2 catalogue. As stated in the introduction to the Hip- 
parcos and Tycho Catalogues (ESA 1997) this frame is prac- 
tically identical to the International Celestial Reference Sys- 
tem, and we demonstrate the level of any residual zeropoint 
errors with respect to the ICRF in Section ^j. Furthermore, 
our proper motion zeropoint is taken from extragalactic ob- 
jects, and is therefore also tied to an inertial reference frame. 
The advantages of such a proper motion reference frame are 
discussed extensively with reference to the Lick Observa- 
tory Northern Proper Motion programme (NPM - Klemola, 
Jones & Hanson 1987) and are clearly demonstrated by re- 
cent results from the corresponding southern hemisphere 
survey, the Yale Observatory/San Juan Southern Proper 
Motion program (SPM - Platais et al. 1998; Mendez et 
al. 1999). 



2 METHODOLOGY 

All photographic material for the SSS survey is measured 
on SuperCOSMOS, a fast high precision microdensitome- 
ter. Aspects of the system enabling highly accurate posi- 
tional measurements are described in Hambly et al. (1998) 
and will not be repeated here, except to state that it was 
demonstrated that the machine has the capability to cen- 
troid isolated, well-exposed stellar images on modern fine- 
grained emulsions at a precision of ~ 0.5/im (or 33mas at 
the plate scale of the sky survey Schmidt plates). 



2.1 Global astrometric plate solutions 

Converting the positions of stars and other objects mea- 
sured on a photographic plate into celestial coordinates is 
a classical problem in photographic astronomy and various 
techniques are available to solve it. All the techniques are 
based on having a grid of reference stars, with known celes- 
tial coordinates, distributed over the plate. A transformation 
is then defined between the measured xy positions and the 
celestial coordinates of the reference stars. It is the details 
of calculating the transformation which may vary between 
different techniques. 

The SuperCOSMOS global astrometric software (here- 
after XYTORADEC) uses an algorithm based on the Star- 



link program ASTROM (Wallace 1998). The code imple- 
menting the algorithm is partly borrowed from ASTROM 
and partly consists of calls to routines in SLALIB (Wal- 
lace 1999), the Starlink Positional Astronomy Library. 
ASTROM (and hence XYTORADEC) use a modified ver- 
sion of the standard, traditional technique for converting 
plate positions to coordinates. In this technique the cata- 
logue coordinates of the reference stars are converted to tan- 
gent plane coordinates, adjusted for geometrical distortion 
(caused by the telescope optics and mechanical deformation 
of the plates during exposure) and a linear six-coefficient 
least squares fit is made between the adjusted standard co- 
ordinates and the measured positions. Several authors have 
described this technique, eg. Konig (1962), Taff (1981) and 
Green (1985). In this method a simple, single, global fit 
is made for all the reference stars imaged on the plate. 
The technique used by ASTROM and XYTORADEC dif- 
fers from the traditional method in that the plate centre is 
allowed to vary (that is, included in the fit). An outline of 
the method is as follows: 

First, the point on the celestial sphere at which the 
telescope was pointing when the plate was exposed is estab- 
lished. Reference star coordinates are converted to apparent 
(that is, observed) coordinates in the desired system. Stan- 
dard (tangent plane) coordinates, £77, are computed for the 
reference stars about the tangent point. The standard co- 
ordinates are then adjusted for the geometrical distortion 
where each image centroid is adjusted radially. The adjust- 
ment consists of separately multiplying each £ and r\ by the 
factor r/tanr where r = (£ 2 +'7 2 ) (see Murray 1983 page 
196). 

After applying the geometric distortion a small system- 
atic distortion remains between the adjusted standard co- 
ordinates and positions measured on Schmidt plates. These 
residuals are the well-known 'swirling patterns' seen in the 
residuals of linear fits made to Schmidt data (eg. Taff et 
al. 1992). A set of corrections at a grid of points over the 
plate are prepared beforehand. A grid size of 1 cm on the 
plate (~ 10 arcmin) was chosen since this provides suffi- 
cient resolution to accurately map out the ~ 30 arcmin- 
scale non-linear distortion expected. Different grids are used 
for plates originating from different Schmidt telescopes and 
in different wavebands. These grids are prepared by aver- 
aging the residuals over numerous plates, using this fitting 
procedure (but of course without these corrections). Fig- 
ure [j] shows four examples of the distortion patterns from 
the three Schmidt telescopes providing data for sky surveys. 
For the UK Schmidt J survey, the data are averaged over 
200 plates containing 193866 individual standard star resid- 
uals. The ESO Schmidt data are from 142 plates/107346 
stars while for the Palomar Schmidt first epoch E and sec- 
ond epoch R surveys the data are from 51 plates/47101 stars 
and 9 plates/20814 stars respectively. In all cases, the refer- 
ence catalogue used was the Tycho-2 (H0g et al. 2000). The 
largest non-linear effects are seen for the Palomar POSS- 
I E plates. This has been seen previously, and is probably 
a result of the design of the original plate holders used for 
that survey (eg. Irwin 1994) . Note that the redesign of these 
plate holders has resulted in much smaller non-linearities 
for the second epoch survey - compare Figures [l](c) and (d). 
For each reference star or programme object the appropriate 
correction is derived by interpolating the grid at the posi- 
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Figure 1. Mechanical deformation (or 'swirl') patterns for (a) UK Schmidt SERC-J/EJ plates, (b) ESO Schmidt R plates, (c) Palomar 
Schmidt first epoch E plates and (d) Palomar Schmidt second epoch R plates. In each case, the length of the vectors is scaled such that 
one tick— mark corresponds to one arcsecond. 



tion of the star or object. Simple bilinear interpolation is 
used (eg. Press et al. 1986) . The corrections are then simply 
added to the standard coordinates. 

The standard coordinates, adjusted for mechanical and 
optical distortion, are then fit to the measured positions. 
For a complete description of XYTORADEC see Daven- 
hall (2000). Section ^ illustrates typical results from the 
above procedure. 

2.1.1 Reference catalogue. 

The SGC survey data originally made publicly available on 
the World Wide Web made use of the ACT reference cata- 
logue (Urban et al. 1998). Since mid-2000, the Tycho-2 cat- 
alogue has been available (H0g et al. 2000), and all SGC as- 
trometry was re-reduced with respect to this catalogue. SSS 
data are now routinely reduced with respect to the Tycho- 
2 reference catalogue. This compilation of 2.5 million stars 
reaches magnitudes V ~ 12.5 (90% complete to V ~ 11-5), 
and has positions at mean observational epoch of ~ J1991.5 
from the Tycho catalogue, but includes proper motions via 
a combination of reprocessed Tycho data and first epoch 
positions from the Astrographic Catalogue and 143 other 
ground-based catalogues (see references in H0g et al.). The 
astrometric standard stars provided by Tycho-2 are heavily 
saturated on sky-limited Schmidt plates, and as we show 
in Section there is a demonstrable magnitude systematic 



between bright and faint stars at the level of up to ~ 0.3 arc- 
sec at field edges. Future projects, both ground-based (for 
example, the UCAC project - see Zacharias et al. 2000) and 
space-based, will provide higher density catalogues reach- 
ing fainter magnitudes and will bridge the magnitude range 
(10.0 < V < 15.0) required to eliminate such systematic 
errors from the plate solutions (e.g. Irwin et al 1998). In 
this way, the full-blown internal precision obtainable from 
centroiding on Schmidt plates will ultimately limit the ex- 
ternal accuracy of the astrometry. We note here that the SSS 
survey database is organised in such a way as to make rere- 
duction of the plate astrometry easy to achieve, and when 
higher density catalogues become available this will indeed 
be done. 

2.1.2 FITS World Co-ordinate System for images 
extracted from the SSS database. 

The availability of global astrometric plate solutions for the 
survey data enables any extracted FITS image (see Paper I) 
to contain FITS World Co-ordinate System (WCS) header 
information. As yet, there is no formally accepted standard 
method and keyword set for implementing the WCS. How- 
ever, Greisen & Calabretta (1999) have proposed a stan- 
dard, and we have followed this in our database access soft- 
ware. The WCS is defined locally for each extracted image 
by fitting a linear plate model between tangent plane pro- 
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jected celestial co-ordinates and xy in pixel units. Large- 
scale, non-linear distortions are not relevant to small scales 
of tens of arcminutes. We simply specify the WCS using 

the RA TAN and DEC — TAN prescription along with the CD 

rotation matrix elements (see Greisen & Calabretta f999 
for more details). This set is compatible with, for example, 
the Starlink software collection (eg. image display utilities 
such as GAIA/SkyCAT, Draper 1999). Note also that pre- 
viously defined keywords sets are also written to the FITS 
headers for backwards compatibility. 



2.1.3 R survey data on film 

At the time of writing, it seems likely that the AAO-R sur- 
vey will be completed on Kodak 4415 'Tech-Pan' film. For 
example, a small number (~ 25) of SGC fields have film orig- 
inals measured for the database. Due to its finer emulsion 
granularity, this medium is superior to IIIaF in all aspects of 
astronomical performance (eg. Parker & Malin 1999; Parker 
et al. 2001) with the exception of large scale astrometric 
stability. The 'swirl' patterns seen for glass plates (eg. Fig- 
ure [l]) do not repeat from film to film, but systematics at 
similar scales and of a similar amplitude are nonetheless 
present. Such fields have individual systematic correction 
via the corresponding J plate data. The procedure is simply 
to map out all plate-to-plate errors as a function of position 
between the R film and its corresponding SERC-J/EJ plate 
in a given field using the J plate to define the 'swirl' correc- 
tion mask. 



2.2 Proper motions 

2.2.1 Choice of plate pairs and pairing criterion 

The plates available for the SSS provide several potential 
plate pair combinations for the determination of proper mo- 
tions. Any (or all) of these combinations could be used for 
this task. For the purposes of producing a uniform proper 
motion survey in the southern hemisphere to the highest 
possible accuracy, we chose the J/R combination for sev- 
eral reasons. Each field has (or soon will have) a good qual- 
ity original photograph on fine-grained emulsion reaching 
J ~ 22.5, R ~ 21.0. The POSS-I 'E' and ESO-R material 
available are glass copies, are between 0.5 and 1.0 magni- 
tudes less deep, and neither survey covers the entire southern 
hemisphere; the POSS-I plates have the additional compli- 
cation of not being on the same system of field centres as the 
ESO/SERC atlas. The SERC-I survey lacks depth, reach- 
ing I ~ 19.0. Moreover, the R survey was envisaged as a 
second epoch survey to the J, and efforts have been made to 
ensure a good time baseline between first and second epoch 
(Morgan e t al. 1992; e.g. see Figure 6 later). As we show 
in Section 2.2.4, we can accurately model any colour effects 



introduced in the astrometry as a result of our choice of J/R 
plates. Of course, future enhancements to the SSS may in- 
clude sophisticated analysis of the multiple epoch data avail- 
able from more than two plates per field if there is sufficient 
user demand (e.g. Paper I). 

Figure ^ shows the number of pairings made between 
the J and R plates for several different Galactic latitudes. 
These curves show the same general trends. Below a pairing 
criterion of 1 arcsec, a dramatic fall in the number of pairings 
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Figure 2. The number of objects paired between two plates 
in fields at various Galactic latitudes: b = —32° (solid circles); 
b = —47° (plus signs); b = —65° (open squares) and b = —87° 
(crosses). 



occurs due to global positional errors (see Section |3|) . For a 
pairing criterion of more than 6 arcsec, a sudden and steady 
rise is seen in the number of pairings at the level of many 
thousands at each step. These must be spurious, since we do 
not expect 10 4 — 10 5 high proper motion stars per Schmidt 
field. In fact, the total number of stars having [i > 0.2 arc- 
sec yr _1 is of order 10 over the entire sky (Luyten 1979). A 
cummulative histogram of the fraction of high proper mo- 
tion stars (ie. those having fj, > 0.2 arcsec yr _1 ) made from 
the NLTT catalogue (Figure ||) shows that 90% of those 
catalogued have motions fi < 0.4 arcsec yr _1 . For a median 
epoch difference of 15 yr (see Figure [|), 0.4 arcsec yr~ x cor- 
responds to a shift of 6 arcsec. Hence, in order to pair 90% 
of 'high' proper motion stars, but at the same time avoid 
large numbers of spurious pairings, the pairing criterion has 
been set at 6 arcsec. However, the database access software 
(see Paper I) allows the user the default option of limiting 
pairings to those within 3 arcsec for most applications when 
completeness of high proper motion objects is not an issue. 
This is to further minimise the number of spurious pairings 
for applications that would be sensitive to such contamina- 
tion - clearly, in Figure |^ even between pairing criteria of 3 
and 6 arcsec there is a rise in pairings over and above that 
expected from high proper motion objects alone. As far as 
pairing motionless objects is concerned, a pairing criterion 
of 3 arcsec is very generous in the light of absolute errors in 
position of a a ,s ~ 0.3 arcsec (Section^). 

2.2.2 Plate-to-plate error mapping 

The determination of accurate proper motions requires that 
systematic positional errors as a function of plate position 
be removed from the displacements of all objects measured 
between the first and second epoch plates (eg. EI95). Such 
errors arise from emulsion shifts on the photographs them- 
selves, and also from systematic positional errors during the 
measurement process. Note that in Hambly et al. (1998) 
we described a procedure for minimising systematic mea- 
surement errors arising from the xy table of SuperCOS- 
MOS; nonetheless we follow EI95 exactly and compute one- 
dimensional error functions (ie. errors in x and y as a func- 
tion of x and y) when mapping the plate-to-plate system- 
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Figure 3. Cummulative histogram of the fraction of the total 
number of NLTT stars having a proper motion less than fi arc- 
sec yr _1 . Approximately 90% of NLTT stars have proper motions 
less than = 0.4 arcsec yr -1 . 

atic errors since small residuals are expected to be present 
at the few tenths of a micron level (for example, errors re- 
sulting from pixel placement within the linear CCD used 
to scan the plate). Other systematic errors possibly present 
in the proper motions include those due to positional er- 
rors as a function of magnitude. For example, it is likely 
that the centroids of heavily saturated images off-axis on 
the photographs are displaced relative to those of fainter 
stars owing to off-centre haloes of scattered light, an effect 
compounded by the limited dynamic range of the measuring 
machine. Note however that such effects cannot be simply 
measured from the data and removed since we expect pop- 
ulations of objects having different mean apparent magni- 
tudes to show systematic differences in mean proper motion 
due to Galactic structure (eg. brighter stars will show sys- 
temic motions with respect to fainter stars since they are, on 
average, nearer to the Sun). We argue that such effects will 
be minimised by using the photographic material specified in 
the previous Section since the first and second epoch plates 
were obtained on the same field centres with the same tele- 
scope using the same guide stars at similar times of the year 
and at similar (optimum) airmass. Other effects due to mis- 
alignment of the scanning beam with the plate along with 
different 'depths' of image within the emulsion are likely to 
be present - eg. EI95. In their work such effects were elimi- 
nated by measuring each plate twice, at orientations within 
the machine of 0° and 180°. Taking the mean position be- 
tween the two scans, these errors were removed. The scale 
of the SSS makes scanning each plate twice impractical and 
so systematic errors resulting from such effects will remain 
in these data. 
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Figure 6. Histogram of epoch differences between the SGC Bj 
and R photographic observations. 



ages have zero mean displacement fixes the proper motion 
zeropoint to the extragalactic frame. Note that the image 
classifier is reliable at the level of > 90% down to Bj ~ 20.5 
(Paper II, Table 9) and so only galaxies brighter than that 
magnitude are used in this procedure. Figure ^ shows an ex- 
ample of the distribution in galaxy displacement as a func- 
tion of magnitude for the SGP field 411, demonstrating the 
effectiveness of this procedure. As an independent and ex- 
ternal check on the extragalactic zeropoint, we also exam- 
ined the displacements of all isolated, good quality images 
paired between the SGC survey database and the Veron- 
Cetty & Veron (1998) QSO catalogue. Figure | shows this 
comparison (> 5000 QSOs over the entire SGC region). Ta- 
ble [l| details the mean displacements and RMS errors for 
the QSO sample. Clearly, there are no significant non-zero 
residuals. Moreover, this test illustrates the likely level of 
errors in the stellar proper motions. Figure ^| shows a his- 
togram of the epoch differences between the J and R survey 
photographs within the SGC survey. For a median epoch 
difference of 15 yr, the final column in Table |l] gives the 
implied mean proper motion precision and RMS zero point 
error, in either co-ordinate, as a function of magnitude. Ob- 
viously, some fields have much shorter time baselines than 
the median value and the proper motions in such fi elds will 
be correspondingly less accurate (see Section 2.2.5). 



2.2.3 Zeropoint of the proper motions 

The zeropoint of the proper motion system is fixed using ob- 
jects classed as galaxies (Paper II describes the image classi- 
fication procedure). As described in EI95, once the plate-to- 
plate error mapper has removed systematic errors between 
the first and second epoch plates, in general the mean stel- 
lar displacement will be zero (since stellar images dominate 
the number counts) while that of galaxies will be non-zero. 
Simply applying a global translation to make the galaxy im- 



2.2.4 Correction of systematic effects due to differential 
colour refraction (DCR) 

As noted in EI95 differential colour refraction (hereafter 
DCR) corrections will be required in order that objects of 
different colour to the average have no systematic error in 
proper motions measured from plates of different bandpass. 
Again, these errors must not be measured from the data 
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Figure 4. Distribution of residual displacement in galaxian images in xy as a function of magnitude for the SGP field 411, after 
plate— to-plate error mapping and zeropoint translation between first and second epoch plates. 
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Figure 5. Distribution of residual displacement (between first and second epoch plates) in QSO images in xy as a function of magnitude 
for the entire SGC survey (objects catalogued in Veron— Cetty & Veron 1998). Median zeropoints and RMS scatter as a function of 
magnitude are tabulated in Table |lj 
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Magnitude 


Zeropoint and RMS (arcsec) 


RMS zeropoint 




rang 


c 


rtA 




error^ (mas yr ) 


( TTI £1 s vr ^ 1 


15.5 < B, 7 


< 16.5 


0.015 ±0.132 


0.066 ±0.225 


2.2 


11.9 


16.5 < Bj 


< 17.5 


0.043 ±0.137 


0.052 ±0.148 


2.2 


9.5 


17.5 < Bj 


< 18.5 


0.007 ±0.127 


0.024 ±0.133 


0.8 


8.7 


18.5 < Bj 


< 19.5 


0.004 ±0.132 


0.019 ±0.133 


0.7 


8.9 


19.5 < Bj 


< 20.5 


0.005 ±0.164 


0.027 ±0.152 


0.9 


10.5 


20.5 < Bj 


< 21.5 


0.002 ±0.242 


0.016 ±0.214 


0.5 


15.2 


21.5 < Bj 


< 22.5 


0.044 ± 0.395 


-0.066 ± 0.435 


2.6 


27.7 



1 Assuming median epoch difference of 15 yr. 

Table 1. Residual systematic zeropoint errors in QSO positions between first and second epoch plates from Figure ||, and impl ied proper 
motion errors (both random and systematic) for an epoch baseline of 15 yr, as a function of magnitude (however see Section 3.2). 



since populations of different mean colour (eg. redder young 
disk stars having relatively high metallicity and bluer old 
halo stars having relatively low metallicity) are expected to 
have different systemic motions due to their different ensem- 
ble kinematics. We modelled DCR effects following exactly 
the prescription given in EI95 with the necessary changes 
made for the different emulsion/filter combinations used 
here, with relevant data being taken from Evans (1989) and 
references therein. In Figure [?] we show our computations 
of the coefficient of refraction J? as a function of synthetic 
(Bj-R) colour for the Bj and R passbands along with the 
polynomial fits and linear extrapolations used. 

The DCR corrections are applied as detailed in Equa- 
tions 1 and 2, Section 5 of EI95. The corrections take the 
form of small changes to the second epoch positions as a 
function of colour and as a function of mean galaxy colour 
(since the galaxies define the zeropoint described above). 
Note that in computing the refraction coefficient we assumed 
nominal atmospheric temperature, pressure, humidity and 
unit airmass. Calculation shows that this is good enough for 
our purposes - for example, recomputing R for airmass \ m 
the range 1.0 < \ < 2.0 in steps of 0.1 shows negligible zero- 
point changes in R as a function of colour while the gradient 
remains constant. 



2.2.5 Proper motion error estimates 

In order that some estimate of probable errors in proper mo- 
tion are available when using SSS data, error values in units 
of mas yr -1 are provided with each proper motion. These 
errors are estimated on a field-by-field basis, and therefore 
take into account the small differences in survey plate qual- 
ity as well as the varying time baselines. The estimates also 
take into account the effects of centroiding precision as a 
function of magnitude, individually for stars and galaxies. 
The measurements are made from the global average disper- 
sions in centroids, as a functio n of m agnitude, once the error 
mapping described in Section 2.2.2 has been completed. Of 



course, the assumption in making these estimates is that the 
intrinsic proper motion is small compared to the errors and 
makes negligible contribution to the overall dispersion. This 
is obviously true for galaxies and faint stars, but is an in- 
creasingly poor assumption for increasingly brighter stars. 
The errors are thus overestimated for bright stars. 
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Figure 7. Computational models of R, the coefficient of refrac- 
tion, as a function of (Bj-R) colour. The symbol types are + 
for giants, filled squares for subgiants and open circles for dwarfs. 
The solid lines are least-squares polynomial fits to the dwarf data 
along with linear extrapolations to provide an estimate of R for 
any colour (cf. Evans & Irwin 1995, Figures 11 and 12). 



3 RESULTS AND DISCUSSION 
3.1 Positions 

3.1.1 Residuals of the fitting standards 

Figure ^ shows histograms of the mean RMS error per star 
per SGC plate in Right Ascension and Declination for the 
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Figure 8. Histograms of the mean RMS residual per star per SGC plate in either co-ordinate for reference catalogue (ie. Tycho-2) stars 
(a) SERC J/EJ, (b) SERC ER/AAO-R (excluding data from 4415 film), (c) SERC-I, (d) ESO-R and (e) POSS-I E. 



different survey plate collections included in the SSS survey. 
These plots demonstrate the ability of the plate solutions to 
predict the positions of the reference catalogue stars on the 
plates. It is encouraging that even for the 51 POSS-I E glass 
copies - which are not of the highest quality compared to 
modern atlas copies, and which have a mean epoch of 1954 
(Minkowski & Abell 1963) - the global plate solutions are 
accurate at the 0.3 arcsec level. 

The modal values from these histograms are in the 
range ~ 0.1 to ~ 0.15 arcsec, depending on survey type and 
excluding the POSS-I E data. The SERC-I survey gives the 
best results, with some mean residuals down at ~ 60 mas. 
This result is probably due to a combination of minimal 
proper motion errors (the mean observational epoch of the 
SERC-I survey is close to that of Hipparcos/Tycho) and 
also the fact that the I plates are the least sensitive of the 
surveys, resulting in less saturation and lower intensity scat- 
tered light haloes for the reference catalogue standards. 

3.1.2 Empirical uncertainty in positions relative to the 
International Celestial Reference Frame (ICRF) 

In a recent paper, Deutsch (1999) compared positions of 
a sample of extragalactic radio sources with those derived 
from several currently available wide-field digitised pho- 
tographic survey databases. The comparison sample came 
from the very long baseline interferometric study of Ma 
et al. (1998), including 212 sources defining the ICRF. As 
in Deutsch (1999), we included all defining, candidate and 
'other' sources from Ma et al., since their positional uncer- 
tainties are orders of magnitude smaller than those expected 
from photographic data. Within the SSS survey object cat- 
alogue database, we matched up sources within 3 arcsec of 



the quoted ICRF position but excluded poor quality images 
(e.g. those near very bright stars, step-wedges or plate la- 
bels), images below Galactic latitude |6| = 30° (since we 
are concerned with deep, sky limited plates) and deblended 
images (whose centroids will be, in general, systematically 
wrong - eg. Beard, MacGillivray & Thanisch 1990). Table | 
shows the results of the comparison, where we have split 
the sample into stars (for 'stars' in this discussion we are 
referring to extragalactic objects classed as stellar in SSS 
catalogues - presumably these are QSOs) and galaxies since 
image centroiding is in general a factor ~ 2 less accurate for 
galaxies than for stars (eg. Irwin 1985). We also provide fig- 
ures for stellar samples having magnitudes more than 2.5 m 
above the respective plate limits where the random errors 
due to emulsion noise are reasonably small and constant 
(eg. Lee & van Altena 1983). We give results equivalent to 
those in Deutsch (1999), for each of the survey types, to 
enable direct comparison: columns 3 and 4 of Table show 
the median zeropoint of the respective samples (described by 
columns 1 and 2) with respect to the ICRF; column 5 gives 
the number of objects in each sample; columns 6 through 
to 14 are grouped in triplets, each set of three showing the 
maximum deviation in a, S and R(= {a 2 + 8 2 } 1/2 ) for the 
best 68%, 95% and 99% of the samples. 



Comparing the numbers in Table El with correspond- 
ing values in Table 1 of Deutsch (1999), the accuracy of 
the SGC survey global astrometry is clearly as good as any 
other database for the same survey material. Table |^ also 
clearly demonstrates that, globally at least, the astrometry 
of point-like sources has no zeropoint errors with respect 
to the ICRF larger than ~ 100 mas for many of the plate 
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Table 2. Empirical uncertainty estimates in |6| > 30° celestial positions relative to the ICRF (after Deutsch 1999). For these purposes, 
'stars' refers to point-like extragalactic objects, presumably QSOs. 



collections. In the worst cases (e.g. POSS-I E, ESO-R) the 
global zeropoint errors are no larger than ~ 200 mas. Note 
that these comparison objects have magnitudes in the range 
15 < R < 19 (e.g. Deutsch 1999 Figure Id) while the ref- 
erence catalogue stars have R < 12. Note also that some 
proportion of the apparently random errors listed in Table ^ 
will manifest themselves as local, systematic position errors 
(see the next Section). 



3.1.3 Local position errors as a function of magnitude 
and plate position 

As stated earlier, it is important to investigate systematic er- 
rors as a function of magnitude and plate position since the 
external reference catalogue standards are highly saturated 
on the survey photographs. Off-centre scattered light haloes 
inevitably 'pull' the centroids of these stars (as measured by 
the scanning machine) in a radial direction from the plate 
centre. The direction of this shift on the celestial sphere will 
be discontinuous as one moves from using positional data in 
one field to that in an adjacent field; the amplitude of this 
discontinuity will also vary with magnitude. The discontinu- 
ities in systematic errors with respect to external data are 
likely to be highest in fainter magnitude ranges for global 
astrometric fits based on bright standards. 

Stone, Pier & Monet (1999) provide the means to mea- 
sure the size of such systematic effects via comparison with 
their CCD drift -scan data in the Sloan Digital Sky Survey 
calibration regions. These regions consist of equatorial strips 
of ~ 7.5° in Right Ascension and ~ 3.5° in Declination with 
typical astrometric accuracy quoted as ~ 35 mas at R=10.0 
degrading to ~ 70 mas at R=18.0. Figures ^ and [Hi] show 
the results of this comparison for the R data in the range 
10 < R < 14 and 14 < R < 18 respectively, for calibration 
region 'A' which falls over survey fields 824 to 826. The ex- 
tent of the individual field data are indicated in the plots, 
where a 'seamless' catalogue has been extracted from the 



database as detailed in Paper I. Data from the > 10000 ob- 
jects in common have been binned up in boxes ~ 20 arcmin 
on a side and smoothed/filtered on a 3 x 3 bin scale for the 
purposes of display in Figures ^(e) and |lo|(e) . 

Once again, we find no large systematic offset between 
the SGC co-ordinates and those of Stone et al. (1999). The 
random errors are also consistent with those from Table El 
and are estimated from the distributions in Figures p|jic| (a) 
and (c) to be in the range 0.10 < o a ^ < 0.16 arcsec. Fig- 
ures bltlOj (e) demonstrate the level of systematic differences 
as a function of magnitude (error vectors are scaled such that 
1° = 1 arcsec). A discontinuity in the vector direction can 
clearly be seen in Figure ^(e) when moving from field 826 
to 825; the level of this discontinuity is several tenths of 
an arcsec. In any application having a tight error budget 
(eg. fibre spectroscopy with very small diameter fibres, say 
0.5 arcsec or smaller) one should be careful when using these 
data across a field boundary. We note, however, that for 2dF 
applications (ie. 2 arcsec diameter fibres over 2.0°, Gray et 
al. 1993) these data have the requisite astrometric precision. 

The residuals seen in Figure ^(e) at X ~ 3.0°, Y ~ 
— 1.5° illustrate how potentially misleading the histograms 
in Figure [] can be for estimating true external errors in the 
celestial co-ordinates. Clearly, there are systematic discrep- 
ancies in SGC survey astrometry on scales of a degree or so 
and at the level of up to ~ 0.3 arcsec that are not apparent 
from the residuals of the individual field astrometric fits. We 
reiterate that Figure ^| merely shows how well the astromet- 
ric fits model how the reference stars appear on the plate to 
the measuring machine, and not necessarily where they are 
on the celestial sphere. 

The standard errors of a a _s ~ 0.1 arcsec in Fig- 
ures |9|jic| (a) and (c) are dominated by those for the most 
numerous, faint stars at R ~ 18.5. Towards the plate limits, 
object centroiding becomes less accurate due to emulsion 
noise (eg. Lee & van Altena 1983) and we expect that the 
errors in positions will increase with magnitude (eg. Hambly 
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Figure 9. Comparison against SDSS CCD calibration region 'A' data (field centre a ~ h 40 m , <5 ~ 0°) from Stone et al. (1999) in the 
magnitude range 10 < R < 14. In (e), the systematic error vectors have been scaled such that 1 arcsec = 1°. Dotted lines indicate the 
field data boundaries. 
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et al. 1998, Figure 8) such that at the plate limits of R ~ 21, 
the random errors will be a a ,s ~ 0.3 arcsec. 

3.2 Proper motions 

Platais et al. (1998) describe the Yale/San Juan SPM pro- 
gramme, from which the first release of data contains proper 
motions of 58,880 objects (defined by a set of input cata- 
logues) around the south Galactic pole (SGP). The typical 
accuracies of their data range from ~ 3 mas yr" 1 at V < 15 
to ~ 8 mas yr _1 at V=18.5. 

Figures ^] and ^ show straight-forward comparisons 
between SPM and SGC proper motions in four magnitude 
ranges: 10 < V < 14, 14 < V < 16, 16 < V < 17.5 and 
17.5 < V < 18.5 (the latter is limited by the SPM data). It is 
immediately apparent that there are magnitude dependent 
systematic errors in the SGC proper motions of brighter 
stars. Table |^ quantifies these zeropoint errors for the SGP 
region and also gives the scatter between the SGC and SPM 
data. Given the level of random error quoted in Platais et 
al. (1998) these values should be dominated by errors in the 
SGC measurements. 



Comparing with Figure g and Table |l| there is a consistent 
picture of no magnitude systematics for m > 17.5. However, 
brighter magnitude ranges in Table |l| fail to indicate true ze- 
ropoint errors. This underestimate is due to the fact that the 
data cover the whole SGC survey region, and systematic ze- 
ropoints of different magnitude and sign manifest themselves 
as an increased random error. On the other hand, because 
the SGC-SPM comparison is done in a restricted declination 
range near the SGP, there is a clearer tendency for system- 
atic zeropoints of the same magnitude and sign. The global 
random errors for SGC proper motions are therefore domi- 
nated by magnitude dependent systematics for V < 17, and 
Table ^ indicates this by measuring the apparent SGC-SPM 
scatter without removal of zeropoint errors. Application of 
a global correction to the proper motions to remove these 
zeropoint errors is clearly not possible since their magnitude 
and sign changes over the SGC region. Once again, centroid- 
ing errors will increase with magnitude as one approaches 
the plate limits such that at J,R ~ 22, 21 the proper motion 
errors will be ~ 50 mas yr~ , but with zeropoint error with 
respect to the extragalactic frame of < 1 mas yr _1 . 



4 CONCLUSION 

We have presented a detailed description of the astrometric 
properties of the SuperCOSMOS Sky Survey (SSS)Q Us- 
ing the first release South Galactic Cap (SGC) data, we 
have demonstrated the level of random errors in positions 
and proper motions by comparison SGC survey with exter- 
nal data; moreover we have demonstrated the level of any 
systematic effects in these astrometric parameters as a func- 
tion of magnitude, and in the case of positions as a function 
of field position. For a summary of these results along with 
similar comparisons for other large-scale survey programmes 

t database available online at http://www-wfau.roe.ac.uk/sss 



and a guide to using the SSS database, the reader is referred 
to Paper I. 
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Magnitude fj, a cos(<5) zp fj,g zp and a cr* cos rg\ <*)x 5 

range and a (mas yr^ 1 ) (mas yr -1 ) (mas yr -1 ) 

10.0 < V < 14.0 10.0 9.5 -12.4 8.2 16.4 18.9 

14.0 < V < 16.0 8.5 8.5 -10.3 7.8 13.8 15.9 

16.0 < V < 17.5 3.9 9.6 -5.4 9.1 10.7 11.3 

17.5 < V < 18.5 0.8 14.3 -1.8 14.1 14.2 14.3 

1 'Effective' global random error including zeropoint (zp) error. 
Table 3. Quantitative comparison between SGC and SPM proper motions (see text). 
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